Topological insulators with perfect vacancy superstructure and possible implications 

for iron chalcogenide superconductors 
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Motivated by the newly-discovered intercalated iron chalcogenide superconductors, we construct 
a single orbital tight-binding model for topological insulators on the square lattice with a perfect 
vacancy superstructure. We find that such lattice structure naturally accommodates a non-vanishing 
geometry phase associated with the next-nearest-neighbor spin-orbit interaction. By calculating the 
bulk band structures and the finite stripe edge states, we show that the topological insulator phases 
can be tuned at certain electron fillings in a wide range of the model parameters. The possible 
implications of these results for the iron deficient compounds (A, Tl) y Fe2- x Se2 have been discussed. 

PACS numbers: 74.70.Xa, 73.43.Cd, 73.20.-r, 03.65.Vf 



Topological insulators (TIs) are a new type of quan- 
tum matter which have become evident in electronic sys- 
tems with appropriate lattice structures and the strong 
spin-orbit interaction (SOI) 0, 0] ■ They differ from con- 
ventional band insulators with the non-vanishing Z 2 - 
odd Kramers states characterized by the gapless counter- 
propagating excitations on the boundary of materials. 
Due to the topological nature, such gapless states are ro- 
bust against weak disorder or interactions. In particular, 
when these edge states are in proximity to s-wave super- 
conductors, the Majorana fermions could be realized due 
to the proximity effect 0, Q , and this promises possible 
applications in spintronics and fault-tolerant topological 
quantum computation [B|, || . 

One of the prototype systems for TIs is the Kane- 
Mele model which describes tight-binding electrons with 
the SOI on a two-dimensional honeycomb lattice such 
as in graphene 0, H|- However, the SOI turns out to 
be too small to open a measurable gap in graphene. 
Zhang and his collaborators predicted the existence of the 
two-dimensional TI in HgTe/CdTe quantum well struc- 
tures 0] and it was soon confirmed experimentally by 
observing the quantized residual conductance caused by 
edge states [lfj|. TIs can be also realized in other lat- 
tice systems such as in the diamond lattice model which 



generalizes the 2D honeycomb lattice [TJ| to three dimen- 
sions. As the number of the lattice structures supporting 
the TI are limited in electronic materials, several schemes 
are proposed for realizing the topological quantum states 
in cold-atom systems by taking advantage of the optical 
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lattice in model engineering 

Meanwhile, three-dimensional TIs have been theoreti- 
cally predicted [12| and experimentally realized in several 
realistic materials such as Bi\ x Sb x , Bi 2 Se3, BizTe^, 
Sb 2 Te 3 [H-EDl, and TlBiSe 2 [I3|. Recently, the known 
topological compound Bi 2 Te 3 has been shown to exhibit 
pressure- induced superconductivity with T c ~ 3 K [22 1. 



The real crystal structures in the relevant materials can 



be viewed as, more or less, the layered honeycomb and 
triangle lattices, or the distorted ones like the diamond 
lattice. It remains a challenge to search for the topologi- 
cal superconductors, or the topological compounds with 
intrinsic bulk superconductivity at the ambient pressure. 

Here we suggest that the TIs, or even the topological 
superconductors with moderately high transition temper- 
atures at the ambient pressure, could be possibly realized 
in materials with the crystal structures similar to the 
ones in the 122-type (I4/m) layered iron chalcogenide 
superconductors (A, Tl) y Fe 2 - x Se 2 (with A being the al- 
kali atoms) [2^, [2J]. One prominent feature of these 
newly-discovered materials is that there are iron vacan- 
cies presumably forming the orthorhombic or tetrago- 
nal Fe- vacancy superstructures [HI]. Among them, the 
tetragonal vacancy superstructure (with the Fe-vacancy 
density 20%) is perfect and most stable, as within each 
FeSe layer all Fe-atoms are three-coordinated, exhibiting 
the maximal symmetry with either left- or right-chirality 
but without breaking the in-plane four-fold rotational 
invariance, see in Fig. [lja). This perfect superstruc- 
ture has been clearly observed by a number of experi- 
ments [H-[28|. 

Our suggestion is mainly based on the observation 
that, though being in the square (or cubic) lattice, the 
spin-orbit-induced geometric phase (Berry phase) can be 
naturally protected by the vacancy superstructure. Such 
lattice can be viewed as a bi-partite lattice with each 
unit cell containing a minimal square as a block [29j |. 
In Fig. [2] the geometric phase for an electron hopping 
from the site j to the next-nearest-neighbor (n.n.n.) site 
i with a nearby vacancy is illustrated. Notice that for 
the intra-block n.n.n. hopping, the associated geomet- 
ric phases coming from the two permissible paths, both 
bridged though the Se-atoms in the iron selenides, cancel 
to each other as in the conventional square lattice. Hence 
the chiral vacancy ordered lattices (FigQJb) or FigOTc)) 
preserve the required symmetry as in the honeycomb or 
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diamond lattices 0,0, [n|. 

In the rest of this paper, we will show that the TI 
phases can indeed emerge by tuning either band filling or 
other parameters in a single-orbital tight-binding model 
defined on the square lattice with the perfect vacancy su- 
perstructure (FigOJa)). The tight-binding model Hamil- 
tonian we considered is 



Ht = 5> 



(1) 



Where, c ia and Ci„ are the creature and annihilation op- 
erators for electrons with spin a at the site i. We also 
assume tij = t\ or t% for the nearest neighbor (n.n.) or 
the n.n.n. hoppings within a block, and t%j = t[ or t' 2 for 
the corresponding hoppings between the two n.n. blocks, 
and no hopping happens among other sites, see in Fig. 

mb). 
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FIG. 1: (a)The perfect vacancy superstructure (clockwise). 
The open dots represent the Fe- vacancies and one of the block 
is marked by a dotted square; the x- and y-axes are two direc- 
tions for the \/5 x v5 bi-partite block lattice, (b) and (c): The 
right-handed (clockwise) and left-handed (anti-clockwise) va- 
cancy orderings, respectively. The solid lines are for hoppings 
and the dashed lines are for SOIs. The arrows indicate the 
permissible paths for the inter-block SOIs defined by Eq.(2). 




FIG. 2: An electron at the site i hops to the inter-block n.n.n. 
site j in the presence of a nearby vacancy (with the right- 
handed chirality as shown in Fig[TJb)). This hopping path 
takes a non-vanishing geometric phase induced by the SOI 
bridged by the Se-atoms (big yellow) located below or above 
the Fe-square lattice (small red). 

Next, we introduce the SOI for the two n.n.n. sites 

H\ = iX ^ c \aV<y<y' ' (&i x Zj) c i°' , (2) 



where a = ((J x ,a y ,a z ) are the Pauli matrices, ej and e*j 
are unit vectors along the two bonds in series when the 
electron transfers from j to i. 

By taking the block as a unit cell and choosing the 
coordinate directions as shown in FigQ] [2{|, we are able 
to rewrite the total Hamiltonian H = H t + H\ in the 
momentum space, 



H = YF$h a (k)F a , 



(3) 



where = (c'j 



Acr> °Ber' °C(7> 



J Dg ), with A, B, C, and D 



being the four sites within a block; h^-Qs) is a 4 x 4 trace- 
less Hermitian matrix with the non-zero matrix elements: 
h n2 =U + (t' 2 + i\)e~ lk *, h n3 =h + {t' 2 - i\)e-* k y, 
h tU = t 2 + t[e- lk y, h n3 = t 2 + t[e tk *, h t24 =h + (f 2 + 
i\)e- lk y, /i T34 = t 1 + (tf 2 -i\)e- ik *; and ^(k) = ftf(-k). 

In the absence of the SOI, the bulk band structures 
were solved in Ref.[l^|, exhibiting the Dirac-like disper- 
sions near some specific points in the momentum space 
at the 1/4 or 3/4 electron fillings. Increasing t 2 will gen- 
erally tune the sub-band gaps, and the mid-band gap 
may open when the inter-block hopping is (t 2 ) sufficiently 
weak. In order to investigate the topological properties of 
the system with the SOI, we fix t\ = t[ = 1 for simplicity, 
and consider t 2 = 0, 0.5, 1.0, corresponding to zero, mod- 
erate, and strong intra-block frustrations, respectively. 
We find that for each case either the sub-band gaps or 
the mid-band gap can be opened by tuning A at certain 
electron fillings. 

The topological nature of the insulating phases can 
be then investigated by studying the edge states of the 
strip shaped lattice. The two edges are cut along the 
x direction with sufficient distance (50 unit cells in the 
our calculations). The energy spectra are illustrated in 
Figs[3][5] For t 2 = 0, the three gaps among the four bands 
open gradually with increasing A, as shown in Figs[3lb,c). 
Meanwhile, various in-gap states appear within the bulk 
gaps. Notice that these states are two-fold degenerate for 
each momentum k x . We find that a pair of in-gap states 
cross within the sub-band gap at the 3/4-filling at the 
time-reversal invariant point k x = ir. The similar feature 
is observed at the 1 / 4-filling at another time- reversal in- 
variant point k x — 0. Meanwhile, two pairs of in-gaps 
states cross within the mid-band gap at the half-filling 
near the points k x = ir/2 and 37r/2, respectively. As we 
will show later, all these in-gap states in Figj3]as well as 
the in-gap states in Figs l4|5l are true edge states. 

Of course, the edge states do not necessarily appear 
when t 2 increases. This is illustrated in the case of 
t 2 = 0.5 in FigQJa)(A = 0.2), where no edge state ap- 
pears at the 3/4-filling though the sub-band gap still 
opens. Notice also that sometimes the edge states can 
appear at certain electron fillings but no corresponding 
bulk gap opens. When this happens, the edge states 
may be coupled to the bulk metallic states upon var- 
ious scatterings. We find that the condition for the 
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FIG. 3: (a)A = 0; (b)A = 0.2; (c)A 
t[ — 1 and ti — t' 2 — 0. 



0.6. In all cases, t\ = 1, 



FIG. 4: (a)A = 0.2; (b)A = 0.5; (c)A = 0.8. In all cases, 
ti = 1, ti = 1 , fa = 0.5 and t' 2 = 0. 



gapless edge states to appear within the bulk gap is 
A = A + t 2 — t 2 > 0, as explicitly verified in Fig|4] and 
FigJS] for moderate and strong frustration cases respec- 
tively. In fact, t' 2 modifies A and effectively enhances the 
bulk gap. When A = t 2 — t' 2 , the Dirac-like dispersion 
appears at the 3/4 filling as indicated in FigQJb) where 
t' 2 = 0, t 2 = A(= 0.5). 

The properties of the in-gap or edge states can be fur- 
ther clarified by calculating the amplitudes of the cor- 
responding quasi-particle wavef unctions. The results for 
the in-gap states (marked by the red and blue points in 
Fig. HJc)) are shown in Fig. [SJ where the spin orien- 
tation for each state is also identified. It is shown that 
each in-gap state is spatially localized near one of the 
two edges, while the overlap of the two wave functions is 
exponentially small. Focusing on one of the edges, the 
two states with opposite spin z-components propagate 
along the opposite directions. Hence when the number 
of such edge state pairs are odd, they are robust again 
disorders or weak interactions preserving the time rever- 
sal symmetry Q . In our model, such non-trivial TI phases 
can emerge at the 1 /4 or 3/4 electron fillings in a wide pa- 
rameter region as indicated in Figs. 02b, c), Fig. |UJc) and 
Fig. [3Jb,c), respectively. The existence of non-trivial TI 
phases in our model is compatible with the general clas- 
sification of TIs for systems preserving the time reversal 
symmetry [3(|. 

Finally, we discuss the possible implications of 



our results for the newly-discovered iron deficient 
(TI, K) y Fe2- x Se2 compounds. These compounds are 
essentially multi-orbital materials with sizable 3d-4p hy- 
bridizations, and the Fe-3d electron correlation may be 
moderate. Recent first principle calculations find that 
for x = 0.4 a block-spin Neel antiferromagnetic (AFM) 
ground state can be stabilized in the perfect vacancy su- 
perstructure and a band gap ~ 400 — 500 meV opens 
at the Fermi energy when y = 0.8, suggesting the Fe 2+ 
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valence in the charge neutral parent compound 
The variation in y just shifts the chemical potential but 
does not change the block-spin ordering, implying that 
the perfect vacancy structure is robust. The neutron 
diffraction experiments suggest that such structure is sta- 
ble below 550 K or around [26j|. Therefore, one of the 
most important ingredients for our model to host TIs, 
i.e., the perfect vacancy superstructure which breaks the 
inversion symmetry by its chirality, could remain upon 
various electron or hole dopings (by tuning y). Another 
important ingredient, i.e., the SOI, is intrinsic in various 
chalcogenide-containing materials as in known TIs 13- 
17] . The inter-block frustration hoppings as well as the 
SOIs in the intercalated iron chalcogenides are enhanced 
to due the structure distortion [3jJ,|32j (the corresponding 
Fe-Se-Fe angle is enlarged), so that the condition X + t' 2 > 
t-x can be satisfied. It is also suggested that the mod- 
erate Coulomb interaction can induce the SOI through 
the mechanism of spontaneous symmetry breaking [33j]. 
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FIG. 5: (a) A = 0.5; (b)A 

tl = t'i = t2 — t<2 = 1. 
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FIG. 6: The wave function amplitudes of the edge states cor- 
responding to the red and blue points indicated in Fig. Sc). 
(a)For two spin up states with different momenta; (b)For two 
spin down states with different momenta. 



We can expect that even for relatively strong Coulomb 
interaction, the TI phases could be transformed to the 
topological AFM/Mott insulator phases [3J|. Lastly, re- 
cent experiments have suggested a phase separation 25] 
or co-existence of the block-spin Neel AFM order and 



the superconductivity [26j in (Tl, K) y Fe2- x Se2- Thus 
it is interesting to see whether the TI phases driven by 
the ordered vacancies could be in proximity to or co- 
exist with some kind of superconducting phases in this 
or other similar compounds. Therefore, this class of ma- 
terials provides a possible new route of searching for the 
topological superconductors. 

We would like to thank Minghu Fang and Zhuan Xu for 
helpful discussions. This work was supported by the NSF 
of China, the NSF of Zhcjiang Province, the 973 Project 
of the MOST and the Fundamental Research Funds for 
the Central Universities of China (No. 2010QNA3026). 



[1 
[2 

[3 

[4 
[5 

[6] 
[7] 

[8 

[9 

[10 
[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 

[22 
[23 
[24 
[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 

[33 

[34 



X.L. Qi and S.C. Zhang, Physics Today 63, 33 (2010). 
M.Z. Hasan and C.L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 

L. Fu and C.L. Kane, Phys. Rev. Lett. 100, 096407 
(2008). 

F. Wilczek, Nature Phys. 5, 614 (2009). 

C. Nayak, S.H. Simon, A. Stern, M. Freedman, and S. 

Das Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

A. Y. Kitaev, Ann. Phys. 303, 2 (2003). 

C.L. Kane and E.J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

C. L. Kane and E.J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

B. A. Bernevig, TL. Hughes and S.C. Zhang, Science 
314, 1757(2006). 

M. Konig, et al, Science 318, 766 (2007). 

L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 

106803 (2007). 

H. Zhang, et al, Nature Phys. 5, 438 (2009). 

D. Hsieh et al, Nature 452, 970 (2008). 
Y. Xia, et al, Nature 5, 398 (2008). 
Y.L. Chen et al, Science 325, 178 (2008). 

D. Hsieh et al, Phys. Rev. Lett. 103, 146401 (2009). 
Kuroda et al, Phys. Rev. Lett. 105, 146801 (2010). 
Wu, Phys. Rev. Lett. 101, 186807 (2008). 
B. Shao et al, Phys. Rev. Lett. 101, 246810 (2008). 
D. Stanescu, et al, Phys. Rev. A 79, 053639 (2009). 
Kargarian and GA. Fiete, Phys. Rev. B 82, 085106 
(2010). 

J.L. Zhang et al, PNAS 108, 24 (2011). 

J. Guo, et al, P hys. Rev. B 82, 180520(R) (2010). 

M. Fang, et al, larXiv:1012.5236H (2010). 

Z. Wang et al, arXiv:1101.2059 (2011). 

W. Bao et a UarXiv:1102.0830l (2011). 

V.Y. Pomja kushin et aZ.,larXiv:1102.338"0l (2011'). 

F. Ye et ffl UarXiv:1102.2882l (2011). 

H. Chen, C. Cao, and J. Dai. larXiv:1102.4074l (2011L 
A. P. Schnyder et al, AIP Conf. Pro c. 1134, 10 (2009). 

C. Cao and J. Dai. larXiv:1102.1344l (2011). 

X.W. Yan, M . Gao, Z.Y. Lu, and T. Xiang, 
larXiv: 1102.22151 (2011). 

S. Raghu, X.L. Qi, C. Honerkamp, S.C. Zhang, Phys. 
Rev. Lett. 100, 156401 (2008). 

DA. Pesin and L. Balents, Nature Phys. 6, 376 (2010). 



